Multipole expansions for energy and momenta carried by gravitational waves 
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We present expressions for the energy, linear momentum and angular momentum carried away 
from an isolated system by gravitational radiation based on a spin-weighted spherical harmonics 
decomposition of the Weyl scalar ^£'4. Wc also show that the expressions derived are equivalent to the 
common expressions obtained when using a framework based on perturbations of a Schwarzschild 
background. The main idea is to collect together all the different expressions in a uniform and 
consistent way. The formulae presented here are directly applicable to the calculation of the radiated 
energy, linear momentum and angular momentum starting from the gravitational waveforms which 
are typically extracted from numerical simulations. 
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I. INTRODUCTION 

Gravitational waves are one of the most important 
predictions of General Relativity. Though such gravita- 
tional radiation has not yet been detected directly, there 
is strong indirect evidence for its existence in the form 
of the now famous binary pulsar PSR 1913-1-16, whose 
change in orbital period over time matches to very high 
accuracy the value predicted by General Relativity as a 
consequence of the emission of gravitational waves P, . 
Moreover, there is every reason to believe that the new 
generation of large interferometric observatories (LIGO, 
VIRGO, GEO 600, TAMA) wiU finally succeed in de- 
tecting gravitational radiation within the next few years. 
Also, gravitational waves are one of the most impor- 
tant physical phenomena associated with the presence 
of strong and dynamic gravitational fields, and as such 
they are of great interest in numerical relativity. 

Gravitational radiation can carry energy and momen- 
tum away from an isolated system, and it also encodes 
important information about the physical properties of 
the system itself. The prediction of the gravitational 
wave signal coming from the inspiral collision of two com- 
pact objects has been one of the main themes in numer- 
ical relativity over the years, as such predictions can be 
used as templates that can significantly improve the pos- 
sibility of detection. Nevertheless, it was not until a cou- 
ple of years ago that it was finally possible to perform 
long-term stable numerical simulations of binary black 
hole spacetimes 0, 0, 0, H, I3l- In this way, numerical 
simulations have now reached a stage where it is finally 
possible to extract important astrophysical information 
from the collision of two black holes. For instance, it has 
now become possible to compute the so-called "kick" of 



'Electronic address: ruizm@nucleare s.unam.mxl 



t Electronic address: malcubi@nucleares.unam. mxl 
Electronic address: nunez@nucleares.unam.mx 
§ Electronic address: jryoji@nucleares.unam.mx| 



the final black hole, i.e. the non-zero final velocity of the 
merged black hole, see e.g. [1, 0| and references therein 
(we have not attempted to give a complete list of refer- 
ences here since this field is so active at the moment that 
such an attempt would almost certainly become obsolete 
within a few weeks). This non-zero final velocity is a con- 
sequence of the fact that in non-symmetric situations the 
emitted gravitational waves carry linear momentum with 
them. Similarly, the spin of the final black hole can also 
be computed by subtracting the angular momentum car- 
ried away by gravitational waves from the initial ADM 
angular momentum of the spacetime. 

Historically, there have been two main approaches to 
the extraction of gravitational wave information from a 
numerical simulation. For a number of years the tradi- 
tional approach has been based on the theory of pertur- 
bations of a Schwarzschil d sp acetime developed originally 
by Regge and Wheeler |l(|, Zerilli [ll[, and a number 
of other authors, and later recast as a gauge invariant 
framework by Moncrief [l3|- In recent years, however, 
it has become increasingly common in numerical relativ- 
ity to describe gravitational waves in terms of the com- 
ponents of the Weyl curvature tensor with respect to a 
frame of null vectors [l3l |. using what is known as the 
Newman-Penrose formalism [ij). In this paper we ob- 
tain in a simple way general expressions for the energy, 
linear momentum and angular momentum carried away 
by gravitational waves using the spin-weighted spherical 
harmonics decomposition of the components of the Weyl 
curvature tensor. We also compare these expressions 
with the more common expressions used in numerical rel- 
ativity which are based on the theory of gauge invariant 
perturbations of Schwarzschild [15| (see e.g. [§, [l^). 
Though most of the main ideas and results presented 
here are known, our aim is to collect all relevant expres- 
sions together using a consistent set of conventions and 
definitions. 

This paper is organized as follows. In Sec. [Til we sum- 
marize the general expression for the energy and mo- 
menta of gravitational waves in the transverse-traceless 
(TT) gauge. In Sec. lIIII we obtain the general expressions 
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for the energy and momenta using the Weyl scalars and 
the spin-weighted spherical harmonics. Later, in Sec. IIVI 
we briefly discuss some of the ideas behind the theory 
of perturbations of a Schwarszchild black hole and we 
obtain again the general expression for energy and mo- 
menta in terms of the gauge invariant "master" func- 
tions (^'cvon, *odd) and (Qcvcn, Qodd)- Wc conclude in 
Section |Vl In addition, in Appendix \X\ we discuss some 
important properties of the spin-weighted spherical har- 
monics. 



II. GRAVITATIONAL WAVES AND THE 
PHYSICS IN THEM ENCODED 

Let us start by considering a small perturbation /i^i^ 
to a flat background metric g^^,, such that the full metric 
is given by 



(2.1) 



where ^ 1. Knowing that the gravitational field 

has only two degrees of freedom, we can choose a gauge 
such that there are only two unknown functions, /i+ and 
/i^, that represent the two possible polarizations of the 
gravitational waves. That is, we work in the so-called 
transverse-traceless (TT) gauge in which imposes the 
conditions hoa — h\ — h'^-' \j = (where the bar denotes 
covariant derivative with respect to the background met- 
ric). The perturbation tensor then has the form 
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(2.2) 



where we have assumed that we have a plane wave prop- 
agating along the z axis. Moreover, we will be consider- 
ing only outgoing waves which means that the functional 
dependence on r and t of the metric perturbations is 
of the form /(r — t)/r, so that asymptotically one has 
drh —dth = —h. Furthermore, for large r the radia- 
tion can always be locally approximated as a plane wave, 
so that angular derivatives can be neglected when com- 
pared with radial derivatives (but one must be careful 
when one deals with quantities that do not involve radial 
derivatives such as the angular momentum, see below). 

Having described the properties of the metric pertur- 
bation hf^i,, we proceed to find the flux of energy and 
momentum carried away by the gravitational waves. The 
most straightforward way is to use the Isaacson stress- 
energy tensor (for details see e.g. [l3l)j which describes 
the energy and momentum associated with the gravita- 
tional waves averaged over a few wavelengths using the 
so-called short wavelength approximation (one must re- 
member that in general relativity there is no local expres- 
sion for the energy and momentum of the gravitational 
field). In the TT gauge, and in a locally inertial frame. 



the Isaacson stress-energy tensor is given by 

T,. = ^ {d.hj/^ dXI) , (2.3) 

where ( ) denotes an average over several wavelengths, 
and where a summation over the repeated indices (i,j) 
is implied. Using the explicit form of the in terms 
of and /ix, the Isaacson stress-energy tensor can be 
rewritten as 
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(2.4) 



or equivalently (here and in what follows z will denote 
the complex conjugate of z) 



T — 
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16 TT 



(2.5) 



with H := — ih^ , and where Re(z) denotes the real 
part of z. 

We can now use equation (|2.5p to find the fiux of energy 
and momentum carried away by the gravitational waves. 
Consider first the flux of energy along the direction i, 
which is given in general by T*^'. In particular, the energy 
flux of the gravitational waves along the radial direction 
will then be given in locally Cartesian coordinates by 



dE 
dtdA 



r 



Or 



16 TT 



Re (d^Hd^'H) 



167r 



Re (OtHdrH) 
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with dA the area element normal to the radial direction. 
Using now the relation drh = —h mentioned above, we 
can rewrite this as 

dE 



1 



dt dA 16 TT 



HH) = 



1 



16 TT 



H 



(2.7) 



If we want the total flux of energy leaving the system at 
a given time we need to integrate over the sphere to find 



dE 

lim 



dt r^oo 16 TT 



\H\^dVL 



(2.8) 



where we have taken dA — dfl, with dfl the standard 
solid angle element, and where the limit of infinite radius 
has been introduced since the Isaacson stress-energy ten- 
sor is only valid in the weak field approximation. Notice 
also that we have dropped the averaging since the in- 
tegral over the sphere is already performing an average 
over space, plus the expression above is usually integrated 
over time to find the total energy radiated which again 
eliminates the need to take a time average. 

Consider next the fiux of linear momentum which cor- 
responds to the spatial components of the stress-energy 
tensor T''^ . The flux of momentum i along the radial 
direction will then be given by 



dP, 
dtdA 



T 



1 



16 TT 



167r 



Re (d^HdrH) 



(2.9) 
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where in the last inequality we have used the fact that 
asymptotically diH ~ (xi/r) drH {i.e. we are ignoring 
angular derivatives), and also the relation between the 
radial and temporal derivatives for outgoing waves. The 
vector I introduced above is the unit radial vector in flat 
space, 



x/r — (sin 6 cos ip, sin 6 sin ip, cos ( 



(2.10) 



The total flux of momentum leaving the system will again 
be given by an integral over the sphere as 

^= hm I^d>k\H\^dn. (2.11) 

dt r^co lOTT J 

Finally, we consider the flux of angular momentum. 
Locally, the flux of the i component of the angular mo- 
mentum along the radial direction should correspond to 
EijkX^T'^^ with Cijk the three-dimensional Levi-Civita an- 
tisymmetric tensor (this is just r x p in three-dimensional 
notation). However, in the case of gravitational waves 
this expression is in fact wrong since the averaging pro- 
cedure that is used to derive the Isaacson stress-energy 
tensor ignores terms that go as and it is precisely 

such terms the ones that contribute to the flux of angular 
momentum. A correct expression for the flux of angular 
momentum due to gravitational waves was first derived 
by DeWitt in 1971, and in the TT gauge has the form 
(see e.g. [isj ) 



1 



ijk 



(^Xj dkhab 



dtdA 32.^ ,^.^..^.b + 2S.,hbk)drh-. (2.12) 

The last expression can be rewritten in more compact 
form if we introduce the angular vectors associated to 
rotations around the three coordinate axis. These vectors 
are Killing fields of the flat metric, and in Cartesian co- 
ordinates have components given by = ei^^xj (where 
represents the k component of the vector ^i). In terms 
of the vectors the flux of angular momentum can now 
be written as 



dJj 
~dt 



lim 

r^oo 327r 



(2.13) 



where £(^-hab is the Lie derivative of hab with respect to 
^i, and where we have again taken drh = —dth for outgo- 
ing waves. The appearance of the Lie derivative is to be 
expected on physical grounds, since for rotational sym- 
metry around a given axis {£^ihab = 0) one should find 
that the corresponding angular momentum flux vanishes. 

In order to write the angular momentum flux in terms 
of H as we have done with the flux of energy and linear 
momentum, one must now carefully consider the action 
of the Lie derivative on the perturbed metric. The easiest 
way to do this is to work in spherical coordinates (r, 9, ip), 
in which case the angular vectors have components 



S,x ~ (0, — sin(^, — COS </?cot ^ 
— (0, COS — sin (^cot 0) 
= (0,0,1) . 



(2.14) 
(2.15) 
(2.16) 



It is clear that the vector corresponds to one of the 
vectors of the coordinate basis, which implies that Lie 
derivatives along it reduce to simple partial derivatives. 
To calculate the Lie derivatives in the other directions it 
is convenient to first introduce the two complex angular 
vectors ^± :— ^j, ± iS^y. Furthermore, we will introduce 
an orthonormal spherical basis (e^, ee, e-ip), and define the 
two unit complex vectors e± := 1 / [eg ^ i e^p) . One 
can then show after some algebra that the Lie derivative 
of e± with respect to ^± is given by 



£ 



(2.17) 



Let us now rewrite the metric perturbation hab in the 
TT gauge in terms of the orthonormal basis as: 

hab = h^ [ee)a[ee)b - [e^p)a{e^)b 

+ iee)a{e.p)b + {e^)aiee)b 

= H{e^)a{e^)b + H{e+)a{e+)b. (2.18) 

We are now in a position to calculate the Lie derivative 
of hab with respect to One finds, 

£i±Kb = (e-)a(e-)6 J±H + {e+)a{e+)b J±H , (2.19) 
where we have defined the operators 
J± ■■= £,±da- is e^'f CSC 9 



ztidg cot 9 du, — i s esc I 



(2.20) 



with s the spin weight of the function on which the op- 
erator is acting: s = —2 for H, and s = +2 for H (see 
Appendix IX)) . The last result implies that 



i£^^hab)dth-'' = 

from which we find 

{£^^hab)dth-'' = 

{£i,hab)dth'''' = 



J±H dtH + J±H dtH 

2Rc|j±i7 ati?| , (2.21) 

2Re ^^Xh dtH^ , (2.22) 

2Re ^^JyH^tH^ . (2.23) 



Collecting results, the flux of angular momentum be- 
comes 

(2.24) 



dJj 

dt 



- hm — Re i JiH dtH dfl 

r^oo 167r / 



with the angular momentum operators Ji deflned as 
1 



J+ + J- 

' sin ip do — cos ip (cot 9 + i s esc 9) 
i 



(J,-J_) 



Jz 



+ COS ip de — sin ip (cot 9 + i s esc 9) 



(2.25) 
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Notice that, except for a factor of —ih, these are just the 
quantum mechanical angular momentum operators with 
the correct spin weight [18]. 

We now have expressions for the radiated energy, lin- 
ear momentum, and angular momentum (more details on 
the derivation of these expressions can be found in [l^). 
However, it turns out that the TT coefficients /i+ and 
are not trivial to obtain from a numerical simulation. 
One of the reasons for this is that during a numerical 
simulation one obtains the full spacetime metric, and not 
the background spacetime plus a separate perturbation. 
One therefore needs to relate these metric perturbations 
to some geometric quantities, preferably scalars, that can 
be obtained directly from the available data. In the fol- 
lowing Section, we will concentrate on one such quantity 
that turns out to be ideal for extracting gravitational 
wave information, namely the Weyl scalar 5*4. 



III. THE WEYL SCALAR *4 

A. Definition of *4 

It is well known that the scalars needed to completely 
characterize a given spacetime can be expressed as pro- 
jections onto a null tetrad of the Weyl tensor. These 
scalars, known as the Weyl scalars, have several interest- 
ing properties. For instance, for a suitable chosen tetrad, 
they can be related directly with the gravitational waves 
at null infinity. We will be interested mainly in the Weyl 
scalar 'i'4, which is associated with outgoing gravitational 
radiation and is defined as 



*4 



(3.1) 



with Cap^ii, the Weyl tensor and where fc^ and ifi^^ are two 
vectors of a null tetrad constructed from the orthonormal 
spherical basis in the following way 



k^" 



1 

V2 
1 

1 

72 
1 



(er + e$f) , 



r / 1 



(3.2) 



where the vectors e^, e{f, e^, and are the usual or- 
thonormal basis induced by the spherical coordinates. 
The Weyl scalars (^fQ, ^'i, ^'2, ^3) are similarly defined 
as different contractions of the Weyl tensor with the null 
tetrad. 

It turns out that the complex quantity H = — ih^ , 
defined in the previous Section, can in fact also be written 
in terms of the Weyl scalar ^'4. In order to see this notice 
first that if we are in vacuum far from the source of the 
gravitational waves the Weyl and Riemann tensors coin- 
cide (the Ricci tensor vanishes). Using now the standard 



expression for the Riemann tensor in the linearized ap- 
proximation one can easily show that, for outgoing plane 
waves in the TT gauge traveling along the r direction, 
^'4 takes the simple form 



(3.3) 



while all the other Weyl scalars vanish. This implies that 
for outgoing gravitational waves we can write 



H 



^-4 dt"dt' 



(3.4) 



B. Radiated energy and momentum 



We can now use equation (j3.4p to rewrite the expres- 
sions for the radiated energy and momentum derived in 
the previous Section directly in term of ^'4 [30|. How- 
ever, before doing this it is convenient to project \E'4 onto 
the sphere and describe its angular dependence in terms 
of the spin-weighted spherical harmonics ^F'^™ (see Ap- 
pendix One can easily show that ^'4 has spin weight 
s = —2, so that we can expand it as 



\I'4 



J2 E ^''™(-2r''"(0,0)) 



/=2 m=-/ 



with A''™ the expansion coefficients given by 



A 



Lm 



*4(-2:P''"'(^,0)) dfi 



(3.5) 



(3.6) 



Notice that since we are expanding over the harmonics 
of spin- weight s = —2, the sum over I starts at Z = 2 
(the ^y'"™ are only defined for I > \s\). To simplify 
notation, from now on we will drop the limits on the 
summations and will always assume that the sum over I 
starts at 2, while the sum over m goes from —I to I. For 
summations involving coefficients with indices I' — I ± 1 
and to' = TO ± 1 one should only remember that the 
corresponding coefficients vanish whenever I' < 2 and 
Ito'I > /. 

In order to calculate the total flux of energy leaving the 
system we can now use (|2.8p together with the relation 
between H and ^4, Eq. (j?^ . to find 



dE 
'dt 



lim 



+ 00 16 TT 



*4di' 



dn 



(3.7) 



Using now the orthogonality of the s^''™, we can rewrite 
the radiated energy as 



dE _ ^.^ r2 

dt r^oo 16 TT 



E 



t 2 
A''™ dt' 



(3.8) 



In a similar way, we can express the linear momentum 
radiated in terms of ^'4 and the A''™ coefficients. Using 
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the expression for the radiated momentum, Eq. (|2.1ip . 
and again Eq. I\3A\\ . we obtain 



— hm (t Ij 

dt r^oo 16 TT 



(3.9) 



Substituting now the multipole expansion of '^4 we find 



l,7n I' ,m' 

X / dt' I i''-™' dt' . 



(3.10) 



In order to calculate the integral over the sphere, notice 
first that the components of the radial unit vector li can 
be expressed in terms of scalar (i.e. spin zero) spherical 
harmonics as 



smf cos(p ~ \/ 



27r 
2^ 



sm f sm ip = I \l — 



cos( 



1,0 



(3.11) 
(3.12) 
(3.13) 



We then see that the fiux of linear momentum involves 
integrals over three spin-weighted spherical harmonics. 
Such integrals are given in terms of the Wigner S-lm 
symbols with ^3 = 1, and are also explicitly given in Ap- 
pendix |X1 

Instead of and Py it turns out to be easier to work 
with the complex quantity P+ :— P^ + iPy. After a 
straightforward, but rather long, calculation one finally 
arrives at the following expressions for the flux of linear 
momentum 



dP+ 
~~dt 



lim 

8 TT 



E 

l.rn 



dt'A' 



dt' (a;,„A'^"+i+6,,_„A'-i^"+i 

L Jl+l,m+l 



(3.14) 



167r 



dt'A 



dt' L^A'^'^ + di^^A'-^^^ 



(3.15) 



where the coefficients {ai^ra, i'i,m, ci,m, di,m) are given by 

(3.16) 



^J{1 -m){l + m+ 1) 



1(1 + 1) 



bl.m — TTT 



1 {l-2){l + 2){l + m){l + m-l) 
2l\~ 



{21 ~ 1)(2Z + 1) 



di,. 



2 m 



1 Ul-2){l + 2){l-m){l + m) 



I 



{21- 1){21 + 1) 



(3.17) 
(3.18) 

(3.19) 



Finally, for the flux of angular momentum we go back 
to Eq. ((2724)) to obtain 



= - hm Rei 

at r^oo IQtt 



xJi\^J J 'ifidt"dt'j diVj . (3.20) 

Expressing ^'4 in terms of its multipole expansion and 
integrating over the sphere we now find 



'dt 



lim 

r^oo 16 TT 
t ft' 



I l,m l'm''^~°° 



dt' 



A'^" dt"dt' 
2^' '"'j, (-a^''") dn 



00 •/ —00 

I'm' J ( 



(3.21) 



where the action of the angular momentum operators Ji 
on the spin-weighted spherical harmonics is given in Ap- 
pendix We again obtain integrals that involve prod- 
ucts of two spin- weighted spherical harmonics which sat- 
isfy the usual orthonormalization relations. One can then 
easily find the following expressions for the angular mo- 
mentum carried by the gravitational waves 



dJx 
~dt' 



dJy 

~dr 



dJ^ 
dt 



A^'"' dt"dt' 



>dt', (3.22) 



'^'"'dt"dt' 



(//,m A^'^"^"*^ ~^ fl,-m 
00 ^ 

lim ^ Re/ y f f A' 

(kra - //,-,„ vl'-"-! ) )dt', (3.23) 

lim Im< m / / 

(3.24) 



A''™ dt'dt" 



A''"" dt' 



6 



with 

fi,m ^y{l-m){l + m+l) 

= \/;(/ + l)-m(m + l) , (3.25) 

and where we use the convention that Im(a + ib) ~ ib^ 
for a and b real. These last expressions have also been 
recently derived following a different route by Lousto and 
Zlochower in [23| • The different factor of 1 /4 between our 
expression and the expressions of [20J| is due to a different 
normalization of the null tetrad used to define \l/4. 

IV. BLACK HOLE PERTURBATION THEORY 

In order to relate the expressions for the radiated en- 
ergy and momentum in terms of the Weyl scalar ^"4 to the 
standard ones in terms of gauge invariant perturbations, 
we will present here a brief discussion of some of the ideas 
behind the theory of perturbations of a Schwarzschild 
black hole (a more detailed discussion can be seen in 
e.g. PEIH). 

A. Multipole expansion 

Consider a metric of the form (|2.ip . with the back- 
ground metric g^^ given by the Schwarzschild metric in 
standard coordinates: 

dx^^dx" = fir) de + -1- dr^ + d^'^ , (4.1) 

with /(r) — (1 — 2Af/r). Because of the spherical symme- 
try of the background, it is convenient to think of the full 
spacetime as the product of a Lorentzian two-dimensional 
manifold associated to the coordinates (i, r), and the 
two-sphere of unit radius S'^ associated to {0,(p): 

ds^ = gAB dx^dx^ + r^flab dx^dx^ , (4.2) 

where V,ab is the metric on 5*^: Qat = diag(l, sin^ 6'). 
Here and in what follows we will use upper case indices 
{A, £?,...) to represent the coordinates in M^, and lower 
case indices (a,6, ...) for coordinates in S^. We will 
also distinguish covariant derivatives in the full space- 
time from covariant derivatives in the submanifolds: 
will represent covariant derivatives in spacetime, while 
Da and Da will denote covariant derivatives in and 
5^ respectively. 

We will now consider an expansion of the metric per- 
turbation h^^ in multipoles using spherical harmonics 
y'''"(0, (f). Such a decomposition naturally separates the 
perturbation into even (or axial) modes and odd (or po- 
lar) modes: even modes are those that transform as (—1)' 
under a parity transformation (6*, (/>) (tt — 6*, tt + 0), 
while odd modes transform instead as (—1)'+^. 

In order to decompose /i^^ one further needs to in- 
troduce the scalar, vector and tensor harmonics (these 



tensorial properties refer only to transformations in the 
unit sphere). Scalar harmonics are the usual functions 
yi,m Vector harmonics, on the other hand, come in two 
different types. The even vector harmonics are simply 
defined as the gradient of the scalar harmonics on the 
sphere 

yj'™ := DaY^'"" , (4.3) 

while the odd vector harmonics are 

X^^" := -e.^Ay'^" - -e^e l^^^^i"''™ , (4.4) 

where eab is the Levi-Civita tensor on the two-sphere 
(£01^ = ~M = ^^^^ = sin 6*). 

Similarly, one can define tensor harmonics of even and 
odd type. Even tensor harmonics can be constructed in 
two ways, either by multiplying the scalar harmonics with 
the angular metric flab, or by taking a second covariant 
derivative of the F'^™. However, it turns out that these 
functions do not form a linearly independent set. Instead 
of the DaDbY'^'"' it is better to use the so-called "Zerilli- 
Mathews" tensor harmonics (see e.g. [2lj) defined as 

Z^f :== DaDbY'^"' + + ^abY''^ . (4.5) 

It is easy to show that this tensor is traceless, which 
implies that the resultant functions are linearly indepen- 
dent. Odd parity tensor harmonics, on the other hand, 
can be constructed in only one way, namely 

X'^ = \[DaX'r + DbXl^"^) . (4.6) 

Having defined the vector and tensor harmonics, the 
perturbed metric is expanded in multipoles, and sepa- 
rated into its even sector given by 

fe) ^H'^^Y''^, (4.7) 

\ / even 

fe) ^H^Yt^, (4.8) 

\ / even 

V / even V / 

notice that for the case I = 1 the tensor Z^'™ vanishes, 
and the odd sector is given by 

fe) , - , (4.10) 

\ / odd 

(<) „ = x'r , (4.11) 

\ / odd 

fe) „ = /^''"^af. (4.12) 
\ / odd 

where the coefficients [B^)^, 77^", X''", G"'™, /i^", /i'-™) 
are in general functions of r and t. 

Notice that, since is a constant, both vector and 
tensor harmonics vanish for I = 0. On the other hand, for 
I = 1 the vector harmonics do not vanish, but the tensor 
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harmonics can still be easily shown to vanish from the ex- 
plicit expressions for yi'™. This means that vector har- 
monics are only non-zero for ^ > 1 , and tensor harmonics 
for I > 2. The scalar mode with I = can be interpreted 
as a variation in the mass of the Schwarzschild spacetime, 
while the even mode with Z = 1 is just gauge and can be 
removed under a suitable transformation [22, [H, [23| ■ On 
the other hand, the odd mode with I = 1 can be inter- 
preted as an infinitesimal angular momentum contribu- 
tion, i.e. a "Kerr" mode (a detailed discussion about 
this point can be found in [221). We will therefore only 
be interested in modes with I > 2, just as it happened 
with the expansion of ^'4. 



B. Gauge invariant perturbations 

As we have seen, even parity perturbations are char- 
acterized by the coefficients (ff^^, i?^™, if''™, G''™). 
These coefficients are clearly coordinate dependent, and 
in particular change under infinitesimal coordinate trans- 
formations of the form x^^^x^^+^f^, with |^^| <^ 1. 
However, it turns out that one can construct gauge in- 
variant combinations of coefficients. Two such invariant 
combinations are (for details see e.g. [HI, [HjlHl) 



Lm 



K' 



Lm 



H 



with e 



I An 



= H' 



AB 



l.m 



DASr 



, ) L.rr 

DbSa 



(4.13) 
(4.14) 



'/2. In terms of the gauge 

invariant perturbations and Hab defines the 

so-called "Zerilli-Moncricf master function" as 



The second equality shows that in fact depends only 
on Ha and not on h (the contributions from h cancel when 
contracted with the e^^), but it is nevertheless gauge in- 
variant. Again, using the perturbed vacuum Einstein 
field equations one can show that obeys a simple 

wavelike equation known as the "Regge- Wheeler equa- 
tion" [13. 



C. Gravitational radiation in the TT gauge 

The advantage of using the gauge invariant quantities 
introduced above is that they can be naturally related to 
gravitational waves in the TT gauge. 

Consider first the even parity perturbations. If asymp- 
totically we approach the TT gauge then we will find that 
Hab and hAb decay much faster than hat- According to 
the multiple expansion we can then ignore the coefficients 
i?^™ and iJ^", and concentrate only on X'^™ and G'^". 
Considering a local orthonormal basis aligned with the 
angular directions one then finds that 



(h 



+ \ l.'rn 
even 



[h 



X \l,m 
/even 



QLn 



G' 



sm 



(4.18) 
(4.19) 



On the other hand, from the traceless condition one finds 
that K = {). In that case the Zerilli-Moncrief func- 
tion (|4.15[) simplifies to ^p™„ = rG''™, which implies 
that the contribution from even perturbations to the TT 
metric functions can be written in terms of 5* even as 



2r 



AB 



vDaK 



Lm 



(4.15) 



where L l{l + 1), A := (/ - 1)(/ + 2) + 6M/r, and 
rA '■= Dat. This quantity can be shown to obey a simple 



wavelike equation known as the "Zerilli equation" [llj, 
though we will not go into such details here. 

One can also construct a gauge invariant quantity for 
the case of odd perturbations in the following way (again, 
see e.g. [Mj, [H, l^) 



" l.m 



1.771 



-t^Da 
2 



(4.16) 



In terms of li^"' one now defines the "Cunningham- 
Price-Moncrief master function" as 



2re 



AB 



odd 



{l-l){l + 2) 

2re^^ 
{I -I) [1 + 2) 



DAhr, 



DaK 



l.jn 



2rA ri 

'T'I 

r 

2rA 



l.m 



(4.17) 



ih+) 



l,m 
even 



even I r7l,m 



2r 



sm 



r 92 1 



Y 



l.m 



(4.20) 



[h 



X \/,m 
/even 



even ( 



sin 9 



sin 6 



d 

— — cot ( 



r''™, (4.21) 



where we used the fact that 90,1" '™ — imY'- 



Consider now the odd perturbations. We find. 



2r2 



Sip 



sin 9 



(4.22) 
(4.23) 



We now need to relate /i''™ to '^^^2d- -'^^ ^^^^ case, how- 
ever, we cannot just ignore /i^™ in favor of /i''™, since 
from the definition of '^l^Zi that it in fact depends 
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only on hj^ and not on /i''*". Nevertheless, in the TT 
gauge these quantities are related to each other. In or- 
der to see this, consider the transverse condition on /i^a, 
V'^/i^a = 0. Using the multipole expansion, and calcu- 
lating explicitly the divergence of X'^l^, one finds that 
this condition implies 



D^'lr^h'r] = i (/-!)(/ + 2) /i'' 



(4.24) 



Remembering now that in the TT gauge one has the extra 
freedom of taking h^i, to be purely spatial, plus the fact 
that asymptotically the metric gAB is just the Minkowski 
metric, the previous expression reduces to 



(4.25) 



In the same limit one can also rewrite expression (|4.17p 

(4.26) 



for ^-^'^ as 



V]/' 



2r 



idd 



(;-!)(; + 2) 

Collecting results we find that 



Since for an outgoing wave we have St/i^'™ ^ 
we can integrate the above expression to find 



l.m 
)dd 



(4.27) 



(4.28) 



We can then rewrite the odd metric perturbations as 



)odd 



odd 




vl.7n 


2r 


sin^e* 




f im \ 


■ d 


r 


\sme J 




odd 






r 


I sin6' ] 





(4.29) 



(id 



inn + n 1(1 



1] 



Y 



(4.30) 



The full TT coefficients and then take the form 
2r \ sin" 6* / 

(4.31) 




(4.32) 



Notice that while the functions (^'evcn' *cvon) 
general complex, the TT coefficients and must 
be real. Using the properties of the spherical harmonics 
under complex conjugation one can easily show that this 
implies 



(-1)™* 



l. — m 
even 



idd 



4)™ V[/'' 

/ Oi 



dd 



(4.33) 



One can also rewrite the expressions above in terms of 
spin- weighted spherical harmonics. In order to do this it 
is in fact more convenient to consider the complex com- 
bination H, for which we find 



H = 



2r E 



'(^ + 2)! 
{1-2)1 



■ dd 



2Y 



l.rn 



(4.34) 



At this point it is important to mention one very com- 
mon convention used in numerical relativity. We start 
by considering a different odd master function originally 
introduced by Moncrief ^12i] : 



2r^h 



2r^ 



1.771 



r 



(4.35) 



Notice that with this definition is clearly scalar and 
gauge invariant. The Moncrief function just defined has 
traditionally been the most common choice to study odd 
perturbations of Schwarzschild, and because of this many 
numerical implementations use this function instead of 
^odd- possible to show that asymptotically, and in 
the TT gauge, Qyi^ reduces to 

M --dt^^dd- 
One finally introduces the following rescaling 



l,m 
even 



dd 





-2)! 


2{l- 


-2)! 




-2)! 



2(Z-2)!^^ 



(4.36) 

(4.37) 
(4.38) 



In terms of Qc™^ and Qj^d' '-'^^ ^'-'^ finds for the complex 
quantity H: 



H 



V2 



?l.ni 
even 



O''™ dt' 

•^odd 



^y''™ . (4.39) 



D. Radiated energy and momentum 

Just as we did for the case of 5*4, we can now express 
the radiated energy, linear momentum and angular mo- 
mentum in terms of the master functions (^'cvom ^odd) 
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and/or (Qj^^™ 



.(j^gjj, v^Q^). There are two ways in which one 
can do these calculations. One approach is to substitute 
directly equations l|4.34p and (|4.39p into the expressions 
for the radiated energy and momentum. It is much easier, 
however, to start from the expressions for the radiated en- 
ergy and momentum in terms of the multipole expansion 
of 4-4, and relate the coefficients v4''™ ^'T,/,™ ,t,^™ 



to (^f 



even' ^odd 



and(Q:,V^„,Qo'd"d)- 

Comparing the multipole expansion for ^'4 with the 

expansions for the metric perturbations, and using the 

fact asymptotically ^'4 = — -ff, we find 




(4.40) 



Using this one can translate expressions in terms of the 



j^i.m (jjj-gctly into expressions in terms of 
and/or (Qi'™„, Q^; 



even 



' ^ oddJ 



Let us consider first the radiated energy. Using equa- 
tions (|3.8p and (I4.40p . one immediately obtains 



dE 



1 

64^ 
1 

32^ 



E 



E 

l.rn 



{l-2)\ 

Q: 



odd 



l,m 
even 



Ql.m 
odd 



(4.41) 



where, in order to derive these expressions one must use 
the fact that, as a consequence of (|4.33|) . 



E 



^ even odd 



even odd 



= 



(4.42) 



For the linear momentum we start from equations 
(P?Ti)) - (P?T5)) . and again use The calculation is 

now considerably longer, and in order to simplify the ex- 
pressions one must make use several times of (|4.33p (or 
its equivalent in terms of the Q's). The final result is 



dP+ 



2^ [ ''"^ ^even Vodd 



fi+l.m+l I "^evcnVcven '^odd ^odd 

IQ^ [ Vevon Vo. 



dt 



odd 



'; + l,m 1 ^;evcn ^cvcn ' odd ^ odd 



, (4.43) 



(4.44) 



careful when comparing with the first reference, as their 
sums over m go only from to I). 

In a similar way we can obtain expressions for the ra- 
diated angular momentum staring from equations p.22p - 
p.24[) . One now finds, after some algebra. 



l,m ^l,m-\-l 



dd 



Ddd 



dJy 

~dt 



dJz 
dt 



Z,m+1 
even 



+^odd^odd 



327r ^ 

I' 



'"- I ^even ^even ^ ^odd ^odd 



where we have defined 



P' 



odd 



O''™ dt' 

^odd 



(4.45) 

(4.46) 
(4.47) 

(4.48) 



and where again the coefficients /;,m are the same as 
before. Notice that the expressions for dJx/dt and dJy/dt 
are manifestly real. For dJz/dt the term inside the sum 
can be easily shown to be purely imaginary, so that the 
final result is also real. 

Equivalent expressions to the set of equations (|3.8p . 
([XTIl) - ([XT5l) and ([X^ - ljXM)) for the energy and mo- 
menta carried away by gravitational waves in terms of 
^-4, or the the set (ICT]) . g331)-glll) and (|051) - (|071) 
in terms of gauge invariant perturbations, were derived 
by Thorne in [l^. One can directly relate them to the 
results presented here by noticing that 



1 



V2r 
\/2r 



{l+2)jl,-m _|_ ■ ( (1+2) gl.-rr 



, (4.49) 



where in Thome's notation (t — r) is the mass mul- 
tipole momenta of the radiation field, S^'"^{t — r) is the 
current multipole momenta, and where and OS''"™ 

denote the Ith time derivative of these quantities. 

As final comment, notice that in order to simplify the 
notation in this Section we have not explicitly introduced 
the limit of infinite radius. Nevertheless, this limit should 
be understood since all the results are valid only in the 
weak field approximation. 



with the coefficients (a/.m, 6;,m, Q.m, rf/.m) the same as be- 
fore. The last expressions are written only in terms of 
(OovoiiJ Qodd)' trivial to rewrite them in terms of 

(^o'von; ^odd)- Thcsc expressions can be easily shown to 
be equivalent to those recently derived by PoUney et al. 
in 0, and by Sopuerta et al. in [l6j . (but one must be 



v. CONCLUSIONS 

In this paper we have considered explicit expressions 
for the energy, linear momentum and angular momen- 
tum radiated by an isolated system in the form of grav- 
itational waves. Starting from a small perturbation /i^i^ 
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a background metric g^,^, and working in the transverse- 
traceless gauge, we have reviewed the standard expres- 
sions for the radiated energy and momentum based on 
the Isaacson stress-energy tensor. Introducing the Weyl 
scalar \['4 and its multipole expansion in terms of spin- 
weighted spherical harmonics, we have computed explicit 
expressions for the radiated energy and momentum in 
terms of the expansion coefficients A''™. Finally, we have 
also presented multipole expansions in terms of the gauge 
invariant perturbations of a Schwarzschild spacetime. In 
particular, the expressions in terms of the expansion of 
the Weyl scalar ^'4 have the advantage of avoiding the 
need of separating out a Schwarzschild "background" 
in standard coordinates from a numerically generated 
spacetime that can be in an arbitrary gauge. Neverthe- 
less, one still needs to find a suitable tetrad to calculate 
the scalar ^4, and as yet there is no standard procedure 
to find it (though some recent progress has been made on 
this issue, see e.g. |26,])- However, as long as one chooses a 
tetrad that approaches the standard outgoing null tetrad 
in flat space for large r the asymptotic expressions will 
be correct. 

Although most of the expressions derived here are 
known, as far as we know they have never appeared to- 
gether in the literature. We have taken great care of 
explaining the origin of all relevant expressions, using 
clear and consistent conventions and definitions. We be- 
lieve that having all these expressions together will be ex- 
tremely useful when calculating radiated energy and mo- 
menta from numerical simulations of astrophysical sys- 
tems. As a final comment, we have available a Mathe- 
matica script to find the energy and momenta carried by 
gravitational waves, and would be happy to provide it to 
interested readers upon request. 
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APPENDIX A: SPIN- WEIGHTED SPHERICAL 
HARMONICS 

In this appendix we will discuss some important prop- 
erties of the spin- weighted spherical harmonics. Spin- 
weighted spherical harmonics where first introduced by 
Newman and Penrose [l^l for the study of gravitational 
radiation, but they can also be used to study solutions 
of the Maxwell equations, the Dirac equation, or in fact 
dynamical equations for fields of arbitrary spin. 

Consider a complex function / on the sphere that 
might correspond to some combination of components of 
a tensorial (or spinorial) object in the orthonormal basis 



(cr, eg, Cip) induced by the spherical coordinates (r, 0, ip). 
We will say that / has spin weight s if, under a rotation 
of the angular basis {eg, e^) by an angle i/), it transforms 



as / 



f. A trivial example is a scalar function 



whose spin weight is clearly zero. A more interesting ex- 
ample corresponds to a three-dimensional vector v with 

components {v^,v^,v'^). Notice that these components 
are different from those in the coordinate basis (which 
is not orthonormal), and are related to them through 
{v^,v^,v'^) = {v^ ,rv^ ,r sm9 v'^) . Define now two unit 
complex vectors as 



e± := (ee T iSip) /V2 
The vector v can then be written as 



(Al) 



(A2) 



where v 



v"^ , :— {v^ ± iv'^)/\/2. By considering a 
rotation of the vectors {eg,e^) by an angle ip it is now 
easy to see that v'^ has spin weight zero, while the spin 
weight of is ±1. 

The spin-weighted spherical harmonics, denoted by 
sY'''"^{d, ip), form a basis for the space of functions with 
definite spin weight s. They can be introduced in a num- 
ber of different ways. One can start by defining the op- 
erators 



sm 



- d, 



5/ := 



smfc 



d^ - scote] f , 



sm( 



if sin" 



- do- 



sin 9 



+ scotd] f , 



(A3) 



(A4) 



where s is the spin weight of /. The spin- weighted spher- 
ical harmonics are then defined for |m| < ^ and Z > |s| in 
terms of the standard spherical harmonics as 



1/2 



(-1) 



{I + sy. 
{i + sy. 



1/2 



s>0, (A5) 



In particular we have oY'-'"' = F''". The above definition 
implies that 



9 isY'^"') 



[(/_s)(^ + s-|-l)]l/^+ly'■™ , (A7) 
[(/ + s)(^-s + l)]l/^_lr'-" . (A8) 



Because of this, 9 and 9 are known as the spin raising 
and spin lowering operators. One also finds that 

99 (s^''") = - [l{l + 1) - sis + 1)] .y''" , (A9) 
m{sY^^"') = ~ [1(1 + 1) - sis - 1)] sY^'"\ (AlO) 
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so the ^y''™ are eigenfunctions of the operators S3 and 
SS, which are generahzations of the Laplace operator on 
the sphere L^: 



1 



smt 



1 



f 

2 a '-'vJ 



(All) 



For a function with zero spin weight we in fact find that 
L2/ = §3/ = 98/. 

One can also find generalizations of the standard an- 
gular momentum operators for the case of non-zero spin 
weight by looking for operators Jz and J± such that (here 
we are ignoring the factor —ih that normally appears in 
quantum mechanics) [isj 



J,,r''™ = im,Y''^, (A12) 
= i[[l^m){l + l±7n)f''^ sY^'""^^. (A13) 



One then finds that such operators must have the form 

Jz ^ d^, (A14) 
J± = e^^^ \ ±ide - cote d^- is esc 6*1 . (A15) 



The operators for the x and y components of the 
angular momentum are then simply obtained from 
J-i- = Jj. ± ijy^ so that we find: 



J, = ( J+ + J_ ) /2 , Jy = -i[J+- J.) /2 , (A16) 



There are several important properties of the spin- 
weighted spherical harmonics that can be obtained di- 
rectly from their definition. In the first place, one can 
show that the complex conjugate of ^F''™ is given by 



Also, the different jF''™ are orthonormal: 



(A17) 



(A18) 



The integral of three spin-weighted spherical harmon- 
ics is also frequently needed, for example in the calcula- 
tion of the momentum flux of gravitational waves, and 
can be expressed in general as 



(0, ^) (0, ^) 33 r'^'™^ (0, ^) dn 



(2^1 -f 1) (2/2 + 1) (2/3 + 1) 



-,1/2 



47r 

h h h 

-Sl -32 -S3 



^1 h h 

mi 7712 'TIS 



(A19) 



The above expression involves the Wigner S-lm symbols, 
which are related to the standard Clebsch-Gordan coef- 
ficients 7711, Z2, '772|j3, 'tis) through 



h h h 

mi 7772 fTT-s 



{Ii,mi,l2,m2\l3,-m3) 



(A20) 



V2/3 + I 

The Clebsch-Gordan coefficients arise from the addi- 
tion of angular momentum in quantum mechanics, and 
correspond to the coefficients of the expansion of an 
eigenstate \L, M) with total angular momentum L and 
projection M, in terms of a basis formed by the product 
of the individual eigenstates |/i, 777i) |Z2,'772). These co- 
efficients have some important symmetries, though these 
symmetries are easier to express in terms of the 3-/777, sym- 
bols. In particular, the i-lm symbols are invariant under 
an even permutation of columns, and pick up a factor 
(^_\y-i+i2+i3 under an odd permutation. Also, chang- 
ing the sign of all three tti's again introduces a factor of 

(_l)'l+'2+/3 

A closed expression for the Clebsch-Gordan coeffi- 
cients was first found by Wigner (see e.g. [28|). This 
expression is somewhat simpler when written in terms of 
the 3-/771 symbols and has the form 



h h h 



= (-1) 



h-mi r 



7711 "T-2 "T-3 

ih +l2- k)l {h +k- hV- ih + 1-3- liY- ih + m^y. ik ~ 7713)!1 



{h +I2 + I3 + 1)! (^1 + miY. (h - Till)! {I2 + 7712)! ih - m2)! 



E 

fc>0 



kl 



ih + h+mi- k)\ {h - 7711 -h k)\ 



ih -h + h- ky. ih - 7713 - ky. ih -h + "13 + ky 



(A21) 



r 



In the above expression the sum runs over all values of 
k for which the arguments inside the factorials are non- 
negative. Also, if the particular combination of {li,mi} 
is such that the arguments of the factorials outside of the 



sum are negative, then the corresponding coefficient van- 
ishes. A more symmetric (though longer) expression that 
is equivalent to (jA2ip was later derived by Racah [29| . 
but we will not write it here. 
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In the general case equation (|A21[) is rather compH- 
cated, but this is not a serious problem as one can find 
tables of the most common coefficients in the litera- 
ture, and even web-based "Clebsch-Gordan calculators" . 
Moreover, in some special cases the coefficients simplify 
considerably. For example, in the case when mi = li, 
TO2 = h and I3 — mj, = li + 12 one finds 



(A22) 



Another particularly interesting case corresponds to tak- 



ing ^3 = 7713 = (i.e. zero total angular momentum in 
quantum mechanics). In that case we find 



h h 
mi m2 



= (/i,mi,;2,m.2|0,0) 



\/2h + 1 



(A23) 



The cases with I3 = I are also interesting as they ap- 
pear in the expression for the linear momentum carried 
by gravitational waves. One finds, in particular 



h h 1 
mi m2 



h h 1 
mi m2 ±1 



+m2,0 



2 mi Si^^i^ 



v/(2/i + 2) (2^1 + 1) (2/1) 



i {h + m\) {h - mi) j [h - [h + m2) 



h {2li + l){2li-l) 



h {2I2 + 1) (2^2 - 1) 



± 5, 



I {h T mi) {h T W2) 
^1 (2^1-1-2) {2I1 + I) 



, X / (^1 T^i) {li±m2) 

+ '^'-'=+^V2^i(2/i + l)(2Zi-l)+ 



' {h T m2) {I2 ± mi) 
2I2 {2I2 + 1) (2^2 - 1) 



(A24) 



(A25) 
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